Exchange interaction parameters and adiabatic spin-wave spectra of ferromagnets: 

A "renormalized magnetic force theorem" 
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The "magnetic force theorem" is frequently used to compute exchange interaction parameters 
and adiabatic spin-wave spectra of ferromagnets. The interest of this approach is that it allows 
to obtain these results from a non-self-consistent calculation of the (single-electron) band energy 
only, which greatly reduces the computational effort. However, as discussed by various authors, 
this approach entail some systematic error. Here, a "renormalized magnetic force theorem" 
allowing to remove this systematic error without increasing significantly the computational effort is 
presented. For systems with one atom per unit cell, it amounts to a simple renormalization of the 
spin-wave spectrum. This renormalization greatly improves the agreement between calculated and 
experimental Curie temperatures of Fe and Ni. 
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The ab initio study of magnetic interactions in mag- 
netic metals and of interlayer exchange coupling in multi- 
layers has been the subject of a large number of publica- 
tions during the past few years. Many these studies rely 
explicitly or implicitly on the use of the "magnetic force 
theorem" (MFT), which allows to perform the calcula- 
tions non-self-consistently , and by taking into account 
only single-particle energies. This is of great practical 
interest for it reduces the computational effort by several 
orders of magnitude. 

The force theorem was first introduced for the case of 
non-magnetic systems Q . Extensions to the case of mag- 
netic systems have been published by Oswald et al. |2|, 
who focused on the case of magnetic impurities embed- 
ded in a non- magnetic host, and by Liechtenstein et al. 
, who addressed the case of exchange interactions and 
spin-wave spectra of ferromagnetic systems. However, 
as mentioned by Liechtenstein et al. in their 1984 paper, 
their MFT yields exact results only in the limit of infinite 
magnon wavelength. At finite magnon wavelength, their 
prescription entails some systematic error, as emphasized 
recently by various authors 0, 0| . 

In the present Letter I present a "renormalized MFT" , 
which corrects the systematic error entailed by the MFT 
of Liechtenstein et al. Q (hereafter called the 11 bare 
MFT" ) . In the case of periodic systems with a single site 
per unit cell, this results in an extremely simple renor- 
malization of the spin- wave spectrum and Curie temper- 
ature. The Curie temperatures calculated with this pro- 
cedure are found to be in excellent agreement with the 
experimental ones, both for Fe and Ni. 

The problem we want to address here is that of cal- 
culating the energy required to create some static trans- 
verse deviation of the magnetic moments in a ferromag- 
net. The physical motivation for this is that, within the 
adiabatic approximation (i.e., provided that the charac- 
teristic time scale for the dynamics of such fluctuations is 
long as compared with the time scale for the dynamics of 



band electrons), it provides a good approximation to the 
energy of the low-lying magnetic excitations and there- 
fore allows to address from an ab initio point of view the 
thermodynamics of itinerant ferromagnets. In practice 
this results into a mapping of the complicated itinerant 
ferromagnetic system onto an effective Heisenberg sys- 
tem having the same low-lying excitations, but much less 
degrees of freedom. Let n(r) and m(r) = m(r)u(r) be, 
respectively, the charge- and spin-density. For transition- 
metal systems, it is usually justified to neglect (as is al- 
most always done) intra-atomic non-collinearity and fluc- 
tuations (which are generally expected to cost a high en- 
ergy). We therefore have u(r) = ur if r belongs to the 
atomic cell JIr around atom R, and we wish to impose 
a prescribed direction to the unit vector ur of each site 
R 0]. 

The proper way of doing this relies on the "constrained 
density functional theory" of Dederichs et al. |7| , which, 
in the present case, amounts to introduce some local ex- 
ternal field B R , perpendicular to the local magnetic mo- 
ment axis ur, playing the role of Lagrange parameters. 
As we wish to constrain two components of the magnetic 
moment for each site, we have two Lagrange parameters 
per site, which are the magnitude and azimuthal angle of 
the local constraining field. 

In (non-relativistic) spin-density functional theory ||, 
the ground state energy Eq is obtained by minimiz- 
ing the Hohenberg-Kohn functional (HKF) £hk[p] = 
Tq [p] + £ ext [p] + £h [p] + <?xc [p] , with respect to the spinor- 
density p = (jio~q + m • <t), where do is the unit spinor 
and cr the vector-spinor whose components are the Pauli 
matrices. The various terms are respectively, the ki- 
netic energy of a non-interacting system having the same 
spinor-density, the potential energy, the Hartree part of 
the Coulomb energy, and the exchange-correlation en- 
ergy. The constrained ground state energy is found 
by introducing the new constrained HKF (CHKF) 
J-Hk[p> B- 1 ] = £HK[p]+£cons[p, B ], where the constraint 
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term is given by £ cons [p, B^] = - £ R B R • J Qr dr m(r), 

and minimizing J^bk[p,B ] with respect to p(r), and 
to the set B- 1 = {B R }. As we want to calculate 
the energy change associated with some infinitesimal 
transverse fluctuation <5u R away from the ferromag- 
netic configuration, we may attempt to use the vari- 
ational properties of the CHKF in order to minimize 
the computational effort. The difficulty lies in the cal- 
culation of the kinetic energy term, which requires the 
knowledge of the effective one-electron spinor-potential 
w e s[p] = V e eao - B e ff • er of energy eigenvalues £i(w e g[p\) 
(labeled in order of increasing energy) and eigenvec- 
tors which yields p as output spinor-density, i.c, 

P(r) = Et" {|V^(r)| 2 rr + [^(r)tr^(r)] • <x} . The ki- 
netic energy is then given by %[p] = Y,it[ £i(™eff[p\) - 
J dr (riV e g — m • B e ff) . This implicit dependence of 
w c ft [p] upon the spinor-density p is the origin of our dif- 
ficulties. In fact, starting from a trial input potential 
and magnetic field and solving the one-electron problem, 
we are only able to obtain the value of the HKF for the 
output spinor-density, which we did not know a priori, 
but not for some spinor-density chosen a priori. This 
problem can be circumvented by using an auxiliary en- 
ergy functional, first introduced by Harris 9] . The Har- 
ris functional (HF) ^Harris [p, B ], in the present context, 
has the same form as the CHKF Tbk[p, B ], with 7o[p] 
replaced by 

JVci r 
rb,B x ]=£s i (t^B ± ])-/dr (nK' fi -m-B' efr ), 

(1) 

where the new effective spinor-potential w' cff [p, B^] = 
ff • a is defined explicitly in terms of p and 

B x by 

, r Dil , , S{£^[p]+£ H \p]+£^c\p]+£co7m\p^\) 

W CS [PV ](r) S W) , 

(2) 

where the functional derivative is taken at the a priori 
prescribed spinor density. This essential difference be- 
tween w' e g[p, B^] and w c « [p] allows us to calculate ex- 
plicitly the value of the HF for a spinor-density p and 
constrains B^ chosen a priori (provided we know some 
suitable approximation of £ KC [p]), which constitutes a 
great advantage with respect to the CHKF. In view of 
the Kohn-Sham theorem 0, the HF and CHKF obvi- 
ously take the same value Eq for the density p* and 
constrains B^* corresponding to the constrained ground 
state (solution of the constrained Kohn-Sham equation). 
In addition, the HF can be shown to be stationary (but 
not necessarily minimal, in contrast to the CHKF) with 
respect to p and B^ in the vicinity of the constrained 
ground state (p*,B ). The properties of the HF have 
been studied by various authors who found that, in fact, 
it often yields a better approximation of the ground state 



energy, in the vicinity of (p* , B- 1 *), than the correspond- 
ing HKF [13. 

Starting from the ferromagnetic state (with u R = u° 
for all sites R), for which we assume the self-consistent 
density po = uqo~q + tooUo • <x to be known, we perform 
some infinitesimal rotations: ur = u° + <5u R . For the 
ferromagnetic state, the constraining field vanishes every- 
where. For the rotated state, we approximate the energy 
by using the HF ^-Harris [p, B ], evaluated for a trial in- 
put density equal to pi n = noao + mo(uo + <5u) • er, and 
for some trial input constrain B^ (to be specified later). 

If we use the local density approximation (LDA) for the 
exchange correlation term (as is almost always done) , the 
only term in our trial evaluation of the HF which varies 
with Su = {Su R } is the band energy (sum of eigenvalues 
Ei), so that the energy associated with the fluctuation Su 
is 



AE(Su) = J2{^«ffiP^,'Bi}) ~ ^KftiPo,0})} 

i=l 

+0 2 (Sn,Sm,SB ± ), 



(3) 



where O v (x, y, . . . ) is a quantity of order p (and higher) 
in x, y, . . . , and Sn = n* — n , Srn = m* — m , SB = 
B- 1 * Bj^, respectively. We are looking for an expansion 
of AE(6u) to second order in i5u of the form 

AE(Su) = ]T ^4rr' Su r ■ <5u R + 4 (Su) , (4) 

R.R' 

with A RR » = - J RR + (5 RR/ E R „ Jrr"), which defines 
the coupling parameters Jrr' . The sum rule 



^ = Arr ' = 



(5) 



r 



R' 



expresses the fact that the total energy remains invariant 
(within the non-relativistic theory) upon a uniform rota- 
tion of the magnetization. The definition J3J| for Ar R 
implies that it is related to the (static) transverse sus- 
ceptibility X by: 2^4 RR = Mr (x _1 ) rr , Mr/. 

The bare MFT of Liechtenstein et al. [3j amounts to 
make the choice B^ = for the trial input constraining 
fields. They then obtain 



AE(6u)= J2 A nwSn Ii -Su w + 2 (5n,Sm,B ± *), (6) 

R.R' 



with i RR / = - Jrr/ + <5rr. ( J^r" Jrr" ) , and 



Jrr< = -Im f de f dr f dr'B xc (r)G T (r, r') 

T J-oo JO R JO„, 



xB xc (t')G 1 (t',t) 



(7) 



Since Sn and Sm are even with respect to <5u, they are 
generally of second order in <5u; however, the constrains 
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B^* are odd with respect to Su, and therefore generally 
of first order in Su, so that 02(<5n, Sm, B^*) = 02(<5u), 
although 2 (5n,(5m) = 4 (5u). As a consequence, the 
parameters Jrr' (hereafter called the bare exchange pa- 
rameters) are not equal to the true exchange parameters 
Jrr< and entail some systematic error. Note, however, 
that this problem does not occur for the case considered 



by Oswald et al. [2j since the constrains vanish in the 
case they consider. 

Clearly, to calculate correctly the true exchange pa- 
rameters Jrr', we need to take the exact constrains B^* 
(which are still unknown yet) as trial input in our esti- 
mate of the HF. If we do so, we obtain after some algebra, 



AE(5u) = J2 



R,R' 



v4rr<<5ur • <5ur/ + (MrJrr' - /TrrOBr* • 5u R , - -xrr/B r * • Br* 



+ 4 (<5u), (8) 



with Mr = dr m(r) is the magnetic moment of atom R, and where the bare transverse susceptibility XRR' an d 
the exchange-correlation response function K RR i are given by 



XRR = - fde f dr f dr' Im [& (r, r')^(r', r)] , 

77 J-oo JQrt Jn„, 



(9) 



K RR , = - I de dr dr' Im[G T (r,r')S xc (r')G i (r',r) + GHr,r')Bxc(r')G T (r',r)] . (10) 

-OO J Or, J Q R / 



r 



The constrains are obtained by expressing the trans- 
verse moments A/r^ur as resulting from the transverse 
exchange-correlation and constraining fields: A/r<5ur = 



Er( 



^rr<5ur' + xrr'B 



R' 



In order to keep the ex- 



pressions compact we introduce matrix notations: A, A, 
K, X, M are the matrices whose (RR')-elements are, re- 
spectively, ^4rr/, Arr-, Krr, XRR: Mr<5rr'- Insert- 
ing the resulting expression of the constrains, B R * = 

J2r> X -1 (M — K) <5ur/, into the above expression 

L J RR' 

of AE(5u), we finally obtain the exact explicit expres- 
sion of the renormalized exchange parameters: 



1 



(M 



R T )x- 1 (r 



K). 



(11) 



One can easily prove the sum rule Mr = X)r' ^rr' i 
which implies that the constrains B R * vanish for a coher- 
ent rotation of all magnetic moment, and that the sum 
rule (J5J is indeed satisfied. The above result, Eq. (JTTJ, 
the "renormalized magnetic force theorem" , constitutes 
the main result of this Letter. Its importance is that 
it corrects the systematic error introduced by the u bare 
magnetic force theorem" , without increasing significantly 
the computational effort. It therefore provides a gen- 
eral method to calculate exactly the exchange parameters 
Jrr< which is several orders of magnitude faster than a 
fully self-consistent calculation. 

The linearized equation of motion of the transverse 
fluctuations is 



d<5uR dAE n v-^ . n / n 

M R — ^ = 2^— xu° =4^4 RR ,<hi R / xu°, (12) 

R' 



dt 



dSu R 



and the spin-wave energies are given, as usual, by 
the eigenvalues of the symmetric matrix HQ. = 

4M- 1 / 2 Aivr 1 / 2 . 

Some better physical insight into the nature of the 
above renormalization of the exchange parameters can 
be obtained if one performs a simple, but yet quite rea- 
sonable, approximation. Let us define 

Ar = f drB xc (r)m(r) = -L^J RR ,, (13) 

which can be seen as some average of the exchange split- 
ting on site R. The second equality, in the above equa- 
tion expresses a sum rule related to the invariance with 
respect to a global spin rotation. If the magnetization, 
within an atomic cell, is sufficiently rigid, i.e., if intra- 
atomic fluctuations of the spin-density have a large en- 
ergy cost (as is usually the case in transition metals) , one 
has approximately: 



K 



4Jrr' 



RR' 



, and xrr' 



8Jrr' 



(14) 



Ar A r Ar, 

Note that the above relations become exact if do not per- 
form the discretization approximation u(r) — ► ur (see 
footnote |6j). One then obtains the following simple and 
transparent renormalization for the exchange interaction 
parameters, 



A 1-4M 



and for the spin-wave matrix 
HQ. = HQ. ( 



i-A~ l nn 



(15) 



(16) 
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where (A) RR/ = A R( 5 RR/ and HQ. = 4 M^AM" 1 / 2 j s 
the bare (unrenormalized) spin- wave matrix. Noting that 
the bare parameters A RR can be expressed in terms of 
the Stoner parameters (l xc ) RR/ = <5 RR ' A R /(2M R ) as 

2A = Ml xc ^1 — Xl xc ^ M, one easily shows that the renor- 

malization (|15|l can be re-expressed as 



2A = MX~ 



M. 



(17) 



The above expression has the familiar form of the 
random-phase approximation (RPA) result for the trans- 
verse susceptibility; it is important to realize, however, 
that in the present context (as shown by the above deriva- 
tion) this result is exact within the LDA (except for 
the discretization approximation u(r) — > u R , which as 
already indicated, can be removed easily). The present 
approach is therefore formally equivalent to approaches 
based on calculations of the transverse susceptibility || , 
however without the need for self-consistent total energy 
calculations. 

For the particular case of periodic lattices, the above 
equations are most conveniently solved in Fourier space. 
For systems with a single atom per unit cell, and using 
the approximation (|14J) . the renormalization of the ex- 
change parameters leads to the simple rescaling of the 
spin- wave spectrum 



ftw(q) = 



1 - to(q)/A 



(18) 



The above result clearly shows that the bare MFT yields 
correct results only in the limit of long wavelength (in 
particular, it yields the correct spin- wave stiffness D), or 
if the spin-wave energies remain much smaller than the 
exchange splitting A. The above result is in agreement 
with the estimate of the error entailed by the bare MFT 
proposed by Grotheer ||: hui < hu < fi£u/(\ — frw/A) 2 . 

The Curie temperature can be calculated by means 
of the random-phase-approximation (RPA) Green's func- 
tion method |ll| . For periodic systems with a single atom 
per unit cell, one has 

1 



6 1 x ^ 



k B T£ PA M N ^ fiw(q) : 



(19) 



so that, by using approximation (|14f) . one therefore ob- 
tains an extremely simple renormalization of the (RPA) 
Curie temperature: 



knT, 



B-lc 



RPA 



1 - 6 



fc B fg PA 
MA 



(20) 



where T^ PA is the Curie temperature obtained from the 
bare exchange parameters. As seen from Table the 
renormalization of the exchange parameters considerably 



improves the agreement between theoretical and experi- 
mental Curie temperatures of Fe and Ni. 

The method discussed here provides a convenient 
and accurate approach to study the exchange interac- 
tions, spin-wave spectra and Curie temperature of com- 
plex systems such as disordered alloys, ultrathin films, 
nanostructures, dilute magnetic semiconductors (e.g., 
Gai-^Mn^As), etc. 
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ful remarks on this work. 
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TABLE I: Curie temperature calculated within the RPA by 
using the bare (Tc PA ) and renormalized (Tq FA ) exchange pa- 
rameters, as compared with the experimental value (T ^ p ) . 
The bare Curie temperature (T"o PA ) is taken from Ref. p4|: 
the renormalized Curie temperature (Tc PA ) i s obtained from 
Eq. 



system 


T£ PA (K) 


T£ PA (K) 


T£ xp (K) 


Fe (bec) 


950 


1057 


1045 


Ni (fee) 


350 


634 


621 - 631 



